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Abstract 

, An explicit solution of the equation for the classical harmonic oscillator with 

^-j. | smooth switching of the frequency has been found . A detailed analysis of a quantum 

harmonic oscillator with such frequency has been done on the base of the method of 
linear invariants. It has been shown that such oscillator possesses cofluctuant states, 
\ different from widely studied Glauber's coherent and "ideal" squeezed states. 

OO ■ 

Q\ ■ Oscillator models are widely used in many branches of physics, such as quantum op- 

tics, atomic, molecular, and solid state physics. Small vibration of dynamic system can 
be describe in terms of harmonic oscillators in both quantum and classical mechanics. To 
include surrounding influences on the vibration, or to simulate the coupling of the vibra- 
tion with other degree of freedom, one can consider time-dependent parameters specifying 
the Hamiltonian of a harmonic oscillator, for example, mass and frequency. Besides, non- 
stationary oscillator models show essentially nonclassical effects, such as squeezing and 
covariance of their quantum fluctuations. Some examples for such phenomena are: mo- 
tion of one ion in Paul trap, which is precisely described by harmonic oscillator, with 
periodically time-dependent frequency |T|, also Berry phase can be achieved when pa- 
rameters of oscillator undergoes a cyclic change [0-0. Agarwal and Kumar have shown 
that a nonstationary oscillator with linear sweep of the restoring force owns nonclassical 
states H. In the present Letter we study the switching of the frequency of a quantum 



nonstationary oscillator by using the method suggested in |TT], |12| . 
The Hamiltonian of a harmonic oscillator is given by 

H = V H —q 2 , (1) 

where the constants m and Q(t) are mass and the frequency of the quantum harmonic 
oscillator. The case M = M(t) can be reduced to the case M = m ( see for example 



TJ, eq. 122 ). 

We recall the method of linear invariants, developed in series of papers [P|-||12|, which 
we apply to one dimensional case: For each quantum system, described by a quadratic 
Hamiltonian, there is a classical two- dimensional isotropic nonstationary harmonic oscil- 
lator with a Lagrangian (classical) 
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and equations of motion 
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As it is shown in ||15|| , these two real equations are equivalent to one complex classical 
equation of nonstationary harmonic oscillator 



e(t) + n 2 (t)e(t) = 0, 



(4) 



where complex function e(t) = ei(t) + ie 2 (t) completely describes the quantum evolution 
of the system in particular with Hamiltonian ([!]). 

Despite of various investigations of nonstationary harmonic oscillator, smooth switch- 
ing of the frequency for finite interval of time, was not presented in the literature. Here we 
study the behaviour of the nonstationary harmonic oscillator with varying frequency and 
constant mass for finite interval of time. We will show that the found classical solution 
for this case completely determines the quantum evolution of the corresponding quantum 
oscillator too. 

Let we consider switching of the frequency Q(t) in the form; 
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where a and u are real constants, with dimension of time and frequency, respectively. 
The shape of the frequency is shown on figure 1. By direct calculation we can check that 
the function 
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is a solution of the subsidiary classical equation @ of two-dimensional harmonic oscillator 
for the same time- intervals as in (|S|). The complex-conjugate of e*(t) is the other linear- 
independent solution of (^). Figure 2 shows the parametric plot of the real and imaginary 
parts of e(t). 

As it is shown in Appendix 1, Wronsky determinant is one classical integral of motion 
for the equation (||). 

D w (t) = e(t) e*(t) - e(t) e*(t) = -2i. (7) 

We will use this invariant beyond to study the properties of quantum integral of motion, 
especially at the calculation of their commutator. 

So we have solved completely the classical problem of nonstationary harmonic oscil- 
lator with switching the frequency (§]). Let we consider the quantum problem of nonsta- 
tionary harmonic oscillator with Hamiltonian (|IJ), when the frequency is switching by the 
same way as classical oscillator (|5|). The quantum system evolves in time according to 
the Schredinger equation 

ih^9>(q,t)=H(t)*(q,t) 1 (8) 
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where ty(q,t) =< q\^,t > is a function of the state in coordinate representation. We 
define creation and annihilation operators ( f2 = Q(0) ) 



2hmVtr 
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(mQ q — ip) 



(9) 



and it is easy to show that [a, a*] = 1, because quantum-mechanical position and mo- 
mentum operators q and p obey the usual commutation relation [q, p] = %h. According 
to the general method |§,[|l2j presented also in [lf| (see furmula (6) therein), the linear 
invariant which corresponds to our particular Hamiltonian has the form; 



A(t) 



'Q e{t) - 



Operator A'(t) is Hermitian conjugate to A(t) and has the form; 



e(t) )a 
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It is also an integral of motion for the quantum system. Using [a, a)] = 1 and Wronsky 
determinant (^) it is easy to show, that these invariants obey commutation relation 



[A(t),A*{t)] 



(12) 



i.e. their commutator does not depends on time, either. Thus the operators A(t) and 
A\t) satisfy boson commutation relation in every moment t. The Hermitian and anty- 
Hermitian parts of A(t) and A\t) are also integrals of motion, so that one can construct 
a new pair integrals of motion 



Qo(t) 



Po(t) 



h 



-(A\t)+A(t)) 



/Ho 



(Im (e) q 



Im i 



m 



(A\t) - A(t)) = Jn (-mKe (e) q + Tie (e) p). 



(13) 
(14) 



By direct calculation using the found expression for e(t) @, and the equality (f|) one can 
check that = %- - |[Q Q , H] = 0, respectively ^ = % - jr[P , H] = 0. The physical 
sense of these two integrals of motion is that the operators Qo(t) and Po(t) are the initial 
operators of coordinate and momentum of the quantum oscillator, scaled with a factor 
V^e(O) ( e(0) = -e*(0) = ±^lm (e(0)) = ^ ); 



6(0) 

Qo(t) = 4(0) 
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and 



(16) 

in Appendix 1, 



P (t) = P (0) = V fi oe(0)p, 

where the initial condition for e(0) and e(0) were used ( see 
respectively ). 

Qo(t) and Po(t) for an arbitrary quantum system are related with the initial operators 
of coordinate and momentum q(0) and p(0) in the time moment t = 0. The evolutions of 
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q and p for the quantum system are compensated in proper way with IZe and Im parts of 
e(t) and e(t) to keep Qo(t) and Po(t) constant. The special choice of operators A(t) and 
A^(t), which are expressed via solutions of classical harmonic oscillator provides this 



property of the quantum integrals of motion in the method of the linear invariants [12 
In particular this is true for the present quantum oscillator whose motion is determined 
in terms of the found solution (|6|). Indeed, the quantum mean value of the integrals of 
motion Qo(t) and Po(t) are constants, as we will convince later. 

Let us determine the evolution of the first and second moments of the operators q and p 
for the particular oscillator with frequency (|Sp. In |TJ| was proved the following theorem; 
the necessary an sufficient condition one quantum systems to preserve an equality in 
the Schredinger uncertainty relation is the states of the system to be eigenstates of the 
operator A = ua + vaJ + w. Here u, v and w are arbitrary complex numbers, with one 
connection; \u\ 2 — \v\ 2 = 1, and the states are called Schredinger Minimum Uncertainty 
States. The invariants fllCf), ( |Tl"D obey this condition ( the proof is based again on the 
expression of the Wronsky determinant ((?]) ). Let us consider the oscillator in such states, 
i.e. \*5/,t >= \SMUS > and z is the corresponding eigenvalue; A\SMUS >= z\SMUS >. 
At this point, there are two ways to solve the problem; to express q and p in terms of A 
and A\ respectively u, v in terms of e(t), e(t) (|B|) and to find the first and second moments 
for this specific Hamiltonian ([[]). The second way is to use the method of linear invariants 
where these moments are obtained in general form, and to establish the connection with 
particular quantum harmonic oscillator with a switching frequency. We chose the second 



way, following the idea of the method for one- dimensional case from |L2 



Presenting A and A^ from [13] in terms of q and p and solving this system about q 
and p we receive their quantum evolutions, expressed by the solutions of the classical 
two-dimensional harmonic oscillator (e = 6\ + ie2) 
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where a(t), b(t) and c(t) are the time-dependent coefficients in the general quadratic 
Hamiltonian |T2|, |T3[ in front of p 2 , pq + qp and q 2 . Taking quantum mean value of q and 
p we obtain the behaviour of their first moments 



< p >-- 
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For our case (|1]) the Hamiltonian's coefficients are a(t) = b(t) = and c(t) = 2 ^ . 
Consequently, the evolutions of quantum mean value of the coordinate < q > and the 
momentum < p > has the forms; 



< q >-- 



<p>- 



! —(e(t)z* + e*(t)z) 



(21) 
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Here z is the corresponding eigenvalue of the eigenstates \SMUS >, in whose states the 
quantum values are obtained. The phase diagram shown on figure 3 presents the evolution 
of the quantum harmonic oscillator with a frequency (E) ( % — Us, m = lg, a = .5s, 
uj = 1Hz, z — 1 + i0.2 ). The oscillator evolves as an ellipse when the frequency is a 
constant ( t < and t > j-). The ellipse with the smaller horizontal axis corresponds 
to the region t < 0. The bold curve presents the region of the switching frequency. The 
point in the first quadrant corresponds to the quantum mean value of the pair invariants 
(<Qo(t) >,< P (t) >), (0) and dH) respectively. 

c pg = c qp ). The quantum 



There are three second independent moments a q , a p and c qp 



-qp 



deviations o 2 ,oi l as expressions of e(t) and e(t) were found for first time in general case 
of the method of linear invariants in [ffOfl ; 



o*Jt)=ha(t)\e(t)\ l 
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The third second moment, cofluctuation c qp , in terms of u(t) and v(t) was found in p4]| , 
but as an expression of the solution e(t) and its first derivative e(t) of the equation @ in 

0; 

4W '*)lV. (25) 



As far as we have e(t) in explicit form (|^) and a(t) = ^— , b(t) = all second moments for 
the quantum harmonic oscillator with Hamiltonian (P and frequency @ are determined 
completely; 
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Using these formulas we have been plotted the evolutions of the quantum mean values 
a q , a p and cofluctuation c qp , presented on figure 3. We can observe that the quantum 
oscillator possess all kind Schredinger minimum uncertainty states \SMUS >; coherent, 
squeezed and cofluctuant states. The term subfluctuant has been suggested by Glauber 
[16j as more accurate than squeezed, but we will use the traditional squeezed for the states 



with a subfluctuation below the vacuum in a coherent state and cofluctuant for the states 
with nonzero cofluctuation. For example the oscillator is in a squeezed- cofluctuant state in 
the region t < 0, as c qp ^ 0. They are ideal squeezed |17|, [L8| only when the cofluctuation 



Cqp = 0. In this situation the Schredinger uncertainty relation devolves in Heisenberg one. 
In the region of the switching frequency the state is squeezed-cofluctuant, too. The same 
situation exist in the region t > except in the time moments 



(n = 1, 2, 3, ...) were the state is coherent ( no squeezing, no cofluctuations - dimensionless 
fluctuations are equal; mfl(t n )a 2 = a 2 /mfl(t n ) = h/2, z is equals to Glauber's complex 
number z = a). The analysis shows that in the moment t = for the initial conditions 

e (0) = \J 1+ 0T anc ^ ^(0) = i\j i+au} ^ ne oscillator is in the ideal squeezed state. 

An other possible description of the quantum problem is based on the important 
Wigner function, which we are going to concern here. The all information on the quantum 
system is contained in the time-dependent density operator p(t), satisfying the conditions 
of hermiticity, normalization and nonnegativity. 

For pure states, the density operator which is projector on the state \^/,t >, i.e. 

p(t) = >< (30) 

satisfies the extra condition 

P 2 (t)=p(t), 

and consequently, 

Trp 2 (t) = 1. (31) 

The density operator ( |30| ) obeys the equation for invariants = ( which differs from 
Heisenberg equation of motion, see for example |y| ) 

S -m + l. m ,H { t)} = (32) 
Instead of density matrix one can consider Wigner function pi . Such formulation of the 



quantum problem by means of functions on phase space is very convenient for quantum 
systems having classical analog. In Weil representation mean value of every physical 
variable is an integral of this variable with the distribution function over the all phase 
space. This procedure is a full analog to the classical one. 

The Wigner function is Fourie-transformation of the coordinate representation of den- 
sity operator. 

roc y y 

W(q,p)=J^ p(q+ -,q~ -)exp(-ipv)dv (33) 

where p(q,q',t) =< q\p(t)\q' > is the density matrix. 

With the help of ([3~0"1), one can calculate the Wigner function (^) in states \SMUS >; 

W(q,p,e(t),^l)= (34) 
2 r 2 r 1 1 

^ exp 7 r <P>f ~ 2ct <ip(p- <P>)(q- <<?>)+ o-J(g- < g >) 2 \ j , 

where for the determinant of the variances matrix det(a(t)) we have used the following 
expression; 

h 2 

det(a(t)) = ay q - c 2 qp = — . (35) 
As far as a p , a q and c qp (|26| - |28D are functions of the solution of the equation (|J), we have 



obtained the evolution of Wigner function in terms of e(t) and ' e }." for our particular case 
of frequency @. On the figure 5 Wigner function in the time moment t = is shown. 
It corresponds to the initial point of the bold curve on figure 3, where the phase-space 
diagram of the quantum harmonic oscillator is presented. 
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In conclusion, an explicit solution of the equation for the classical harmonic oscillator, 
with smooth switching frequency has been found. A detailed analysis of a quantum 
harmonic oscillator with such frequency has been done on the base of the method of 
linear invariants. It has been shown that such oscillator also possesses cofluctuant states, 
different from widely studied Glauber coherent and "ideal" squeezed states. It has been 
found the evolution of the Wigner function of such quantum oscillator with a switching 
frequency. 

The author is grateful to D.Trifonov for valuable discussions. This work was partly 
supported by Bulgarian Scientific Foundation, grant number F-559. 

1 Appendix 1 

Here we will calculate the Wronsky determinant for the classical equation of the two- 
dimensional oscillator (|4]). To derive the Wronskian we need to know the first derivative 
of the function e(t) in the three regions; 



de(t) 
dt 



V l+auj \ Y (1+aw) 1 V l+auj \ V (l+aw) 1 



JO _L +a(cos(ljt)) ^ / n J.\ I *\ 

—a io sin (2a; t) + %) 



v / ^+ Q ( cos (^*)) 2 



e'^ \ — yuj(l — au) sin(a;\/l — auj (t — ^j)) + i-Jui cos(to , v / l — a ^ {t — ^j)) 

(36) 

By direct calculation using @ and (|36| ) we receive the Wronsky determinant; 



D w (t) 



e(t) e*(t) 

de(t) de*(t) 
dt dt 



e{t)e*{t)-e*{t)e{t) = -2i. (37) 



Obviously Dw(e(t), e(i)) is an ( classical ) integral of motion for equation (|). 
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